Monte Carlo study of multicomponent 
monolayer adsorption on square lattices 



G. D. Garcia F. O. Sanchez- Varretti*^'^, F. Bulnes 
A. J. Ramirez-Pastor ^'-"^ 



a 

5 



^Dpto. de Fisica, Instituto de Ftsica Aplicada, Universidad Nacional de San Luis - 
CONICET, Chacabuco 917, 5700 San Luis, Argentina. 

^ Universidad Tecnologica Nacional, Regional San Rafael, Gral. Urquiza 314, 5600, 

San Rafael, Mendoza, Argentina. 



Abstract 

The monolayer adsorption process of interacting binary mixtures of species A and 
B on square lattices is studied through grand canonical Monte Carlo simulation in 
the framework of the lattice-gas model. Four different energies have been considered 
in the adsorption process: 1) eo, interaction energy between a particle (type A or 
B) and a lattice site; 2) waa, interaction energy between two nearest-neighbor A 
particles; 3) wbb, interaction energy between two nearest-neighbor B particles; 4) 
WAB = wba, interaction energy between two nearest-neighbors being one of type A 
and the other of type B. The adsorption process has been monitored through total 
and partial isotherms and differential heats of adsorption corresponding to both 
species of the mixture. Our main interest is in the repulsive lateral interactions, 
where a variety of structural orderings arise in the adlayer, depending on the in- 
teraction parameters {waa, wbb and was)- At the end of this work, we determine 
the phase diagram characterizing the phase transitions occurring in the system. A 
nontrivial interdependence between the partial surface coverage of both species is 
observed. 
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1 Introduction 



The adsorption process of mixture gases on solid surfaces is a topic of great 
interest not only from an intrinsic but also from a technological point of 
view, due to its importance for new developments in fields like gas separa- 
tion and purification p!|2|3|l] . Although this problem has been theoretically 
[5|6|7f8|9|10|ll|12] and experimentally [TT|12|13|14|15] studied for many years, 
some aspects are still unclear being necessary to reach a better understand- 
ing about the behavior of the adsorbate during the adsorption process of the 
mixture. 



As in any adsorption process, a complete analysis of the behavior of gas 
molecules under the influence of an adsorbent requires the knowledge of the 
forces of molecular interactions [T|16fl7fl8] . In other words, the description 
of real multicomponent adsorption requires to take into account the effect of 
the lateral interactions between each species of the mixture. An exact treat- 
ment of this problem, including ad-ad interactions, is unfortunately not yet 
available and, therefore, the theoretical description of adsorption relies on sim- 
plified models [1]. One way of overcoming this complication is to use Monte 
Carlo (MC) simulation method [T9|20f21f22f23] . MC technique is a valuable 
tool for studying surface molecular processes, which has been extensively used 
to simulate many surface phenomena including adsorption [23], diffusion [25] , 
reactions, phase transitions [25], etc. 



In this line of work, a previous article was devoted to the study of the adsorp- 
tion of interacting binary mixtures on triangular lattices [27]. In Ref. [|27]]. 
Rinaldi et al. obtained adsorption isotherms and differential heats of adsorp- 
tion corresponding to both species of the mixture, for different values of the 
lateral interactions between the adsorbed species. An unusual feature was ob- 
served when (i) the lateral interaction between A and B particles was different 
from zero, and (ii) the initial concentration of B particles was in the range 
[0.3,0.5]. In these conditions, the A particles adsorbing on the lattice expel 
the B adsorbed particles; then, the partial A [B] coverage increases [decreases] . 
During this regime the number of desorbed particles is greater than the num- 
ber of adsorbed particles which results in a decreasing of the total coverage 
that occurs for a wide range in values of chemical potential where the slope 
of the adsorption isotherm is negative. The behavior of the system was fully 
explained through the analysis of the phase diagrams for order-disorder tran- 
sitions occurring in the adsorbed layer. 



Because the structure of lattice space plays a fundamental role in determining 
the statistics of mixtures, it is of interest and of value to inquire how a specific 
lattice structure influences the main thermodynamic properties of adsorbed 
mixtures. In this context, the objectives of the present paper are (1) to extend 
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the previous work to square lattices using the same techniques developed in 
Ref. [[27]] and (2) to study the effect of the lattice structure on the adsorp- 
tion of interacting binary mixtures. For this purpose, multicomponent gases 
adsorbed on square lattices are studied by using grand canonical ensemble 
MC simulation. The process was analyzed by following total and partial ad- 
sorption isotherms as well as differential heats of adsorption corresponding to 
both species of the mixture. The detailed behavior of these quantities will be 
shown to be directly related to the phase diagrams of the system. 

The outline of the paper is as follows. In Sec. II we describe the lattice-gas 
model and the simulation scheme. In Sec. Ill we present the MC results. 
Finally, the general conclusions are given in Sec. VI. 



2 Lattice-Gas Model and Monte Carlo Simulations 



The adsorptive surface is represented by a two-dimensional square lattice of 
M = L X L adsorption sites, with periodic boundary conditions. The sub- 
strate is exposed at a temperature T to an ideal gas phase consisting of a 
binary mixture of particles A and B with chemical potentials fiA and /i^, re- 
spectively. Particles can be adsorbed on the lattice with the restriction of at 
most one adsorbed particle per site and we consider a nearest- neighbor (NN) 
interaction energy wxy (X, Y = A, B) among them. The adsorbed phase is 
then characterized by the Hamiltonian: 



1 

^ = o H H [wAASci,ci,l + WBBSci,ci,-l + Wab + + 

i le{NN,i} 
AI M 

+eo + - + f^BSc„-i) (1) 

i i 

where Cj is the occupation number of site i (cj = if empty; Cj = 1 if occupied 
by A and q = —1 if occupied by B); / € {NN,j} runs on the four A^A^ sites 
of site i] the (5's are Kronecker delta functions and eo is the interaction energy 
between a monomer (type A or B) and a lattice site. In this contribution, the 
chemical potential of one of the components is fixed throughout the process 
{Hb = 0), while the other one (/i^i) is variable, as it is usually assumed in 
studies of adsorption of gas mixtures In the actual implementation of the 
model eo was set equal to zero, without loss of any generality. 

With respect to the computational simulations, the Monte Carlo procedure 
used has been discussed in detail in Ref. [[27]] and need not be repeated here. 
In this case, the first 10^ Monte Carlo steps (MCS) were discarded to allow 
equilibrium, while the next 10^ MCS were used to compute averages. 
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3 Results and discussion 



The computational simulations have been developed for square Lx L lattices, 
with L = 96, and periodic boundary conditions. With this lattice size we 
verified that finite-size effects are negligible. Note, however, that the linear 
dimension L has to be properly chosen such that the adlayer structure is not 
perturbed. 

In order to understand the basic phenomenology, we consider in the first place 
the case of single-gas adsorption. This was achieved by making f-iB ^ — oo. 
Fig. 1 shows the behavior of the adsorption isotherms and the differential heats 
of adsorption for different strengths of repulsive interparticle interactions. As 
expected, we obtain the well-known Langmuir isotherm passing through the 
point {nA/ksT = 0, 9a = 1/2) when WAA/ksT = (being kB the Boltz- 
mann constant). Two features, which are useful for the analysis of mixed-gas 
adsorption, are worthy of comment: (a) as the NN repulsive interaction is in- 
creased, the coverage at zero chemical potential decreases and asymptotically 
approaches 6a = 0.226; (b) as the repulsive NN interaction passes a critical 
value Wc/kBT ^ 1.763 [2320], a plateau develops in the isotherm at 6'^^ = 1/2 
indicating the appearance of c(2 x 2) ordered phase on the surface. In what 
follows, we consider mixed-gas adsorption but keeping species B at a fixed 
value of the chemical potential /ib/ZcbT = 0. In addition, we have considered 
/c^T = 1 for simphcity, without any lost of generality. 

We start with the case of a binary mixture in the presence of repulsive lateral 
interactions between the particles. The effect of AA interactions is depicted 
in Fig. 2, where WAA/kBT ^ 0, WAB/kBT = and WBB/kBT = 0. We have 
plotted the partial (a) and total (b) adsorption isotherms, and the differential 
heats of adsorption corresponding to the species A (c) and B (d). It can 
be observed in Fig. 2 (a) and (b), that the initial coverage takes the same 
value 6 = 0.5 for all isotherms; this behavior can be explained as follows: 
for fiA/kBT — )■ — oo the A particle coverage is zero while the B particles are 
randomly distributed on the lattice with 6b given by the Langmuir isotherm 
6b = exp{fiB / kBT) / [1 + exp{nA/ kBT) +exp{nB / ksT)], which for ixs/kBT = 
is ^ = = 1/2. 

As discussed above, for a lattice gas of interacting monomers adsorbed on a 
square surface, there exists a critical interaction {wc/ksT ^ 1.763) correspond- 
ing to a phase transition in the adsorbate. The nearest-neighbor coupling, w, 
determines the character of the phase transition: (i) if < (attractive case) 
the system exhibits a first-order phase transition, and (ii) for w > (repul- 
sive case) a continuous order-disorder phase transition occurs in the adsorbate, 
which is observed as a clear plateau in the adsorption isotherms. Following this 
behavior, a well-defined and pronounced step appears in the partial isotherms 
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as the interaction WAA/ksT is increased. Thus, the interaction between A 
molecules determines a c(2 x 2) ordered phase for such particles. Therefore, 
the A isotherm presents a plateau at half coverage [see Fig. 2(a)]. At equilib- 
rium, the B particles occupy half of the empty sites, and the corresponding 
B isotherm presents a plateau at 6b = 0.25; this behavior is a consequence 
of the excluded volume but is not due to the interactions. The total isotherm 
[Fig. 2(b)] is the sum of the partial isotherms, then the plateau appears at 
6 = 0.75. 

In Fig. 2(c), the differential heat of adsorption corresponding to the A 
species is plotted versus 6a- The behavior of the curves can be explained by 
analyzing two different adsorption regimes: (i) for 0.5 < 6 < 0.75 (0 < < 
0.5), the ad-molecules avoid A^A^ occupancy which produces ^ (as for 
waa/^bT = 0), and (ii) for 0.75 < 6 < 1 (0.5 < 6a < 1), the adsorption of 
one more molecule involves an increment cwaa in the energy of the system, 
where c is the lattice connectivity (in this case, c = 4). The maximum in qA 
for 6 — 0.75~ corresponds to the critical coverage at which a dramatic change 
of order takes place in the system (the system passes from the disordered to 
the ordered phase [21] )• A similar situation occurs for the minimum in qA at 
6 0.75+. 

The adsorption of the A species induces an interesting behavior in the B 
isotherm, which also exhibits well-defined steps although the B particles do 
not interact neither with B particles nor with A particles {wAs/kBT = and 
wbe/^bT = 0). This behavior is a consequence of the excluded volume but is 
not due to the interactions. Then, as it is expected, qB is strictly zero over all 
the coverage range [see Fig. 2 d)]. 

We now continue the study of the effect of AA interactions with WAs/kBT = 
and WBs/kBT = 2 (Fig. 3), therefore introducing BB interactions; WAs/kBT = 
2 and WBB/kBT = (Fig. 4) introducing interspecies interactions (and remov- 
ing BB interactions) and, finally, we analyze WAs/kBT = 2 and WBs/kBT = 2 
(Fig. 5) where all interactions are present. 

In the case of Fig. 3, the behavior of the partial [Fig. 3 a)] and total [Fig. 3 
b)] isotherms and the differential heat of adsorption qA [Fig. 3 c)] is similar 
to those in Fig. 2. The main difference is associated to the value of 6b at low 
pressures (which is lower than 0.5), and to the behavior of qB that is not zero 
any more. 

Clearly, for WAA/kBT = 0, A particles are distributed at random and B par- 
ticles start from a low coverage (close to 0.226, as explained in Fig. 1) which 
rapidly decreases as more A particles are adsorbed. Therefore as the total 
coverage increases, B particles interact with each other less frequently and qs 
increases steadily. However, as WAA/kBT becomes sufficiently high so that an 
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ordered phase is formed, a sudden increase in is produced due to a sudden 
increase in the screening effect between B particles produced by A particles. 

An unusual feature is observed in the case of Fig. 4: as /i^/Zc^T increases 
and A particles start to adsorb, B particles are displaced from the surface so 
that the total A + B coverage decreases and shows a local minimum [Fig. 
4 b)]. This effect, which has been previously called mixture effect [27], can 
be explained as follows: as the coverage of A particles is sufficiently high so 
that AB interactions occur, the repulsive character of wab leads to more B 
particles being displaced from the surface than A particles being adsorbed on 
the surface. 



The mixture effect is clearly reflected in the behavior of [Fig. 4 c)] and qb 
[Fig. 4 d)]. In fact, for < < 1/2, A particles do not interact with other A 
or B molecules and, consequently, Qa = Qb = ^ (see insets). The c(2 x 2) phase 
of A particles starts to develop and is completed at = 1/2. For 1/2 < 9a, A 
particles fill the vacancies. In this regime, the coverage of B particles tends to 
zero and does not perturb significantly the adsorption of the A species. The 
important fluctuations in qb are clear signals that the number of B particles 
is practically zero. 

Our simulations show how the competition between two species in presence 
of repulsive mutual interactions reinforces the displacement of one species 
by the other and leads to the presence of the mixture effect. To complete this 
analysis, it is interesting to note that the mixture effect [also called adsorption 
preference reversal (APR) phenomenon] has also been observed for methane- 
ethane mixtures [71l8]|9][T0] . The rigorous results presented in Ref. [[10]] showed 
that, in the case of methane-ethane mixtures, the APR is not a consequence 
of the existence of repulsive interactions between the ad-species, but it is the 
result of the difference of size (or number of occupied sites) between methane 
and ethane. A similar scenario has been observed for different mixtures of 
linear hydrocarbons in silicalite [12], carbon nanotube bundles [33] and metal- 
organic frameworks [3i 



In the case of Fig. 5, the presence of repulsive BB interactions results in a 
initial coverage of B particles close to 0.226. This small fraction of B molecules 
does not perturb significantly the adsorption of the A species and the mixture 
effect disappears. This finding indicates that, in addition to the requirements of 
repulsive lateral interactions between the two species, the amount of particles 
on the surface is essential for the existence of the mixture effect. The rest of 
the figure can be understood following the arguments given above. 

In the following we study the effects of variable AB interactions as shown in 
Figs. 6 to 9. We start with the case where WAA/kBT = WBB/kBT = (Fig. 
6). Here no ordered phases are formed and for sufficiently high wab/^bT an 
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important mixture effect appears. Let us choose for our analysis the curves 
corresponding to wab/^bT = 1. B coverage is initially 0.5. As A molecules 
are adsorbed B molecules are eliminated from the surface in such a way that 
6 decreases. At the same time qa starts at a low value and increases rapidly 
(with a bivaluated behavior) tending to as ^ — ?■ 1, while qs starts near and 
tends to 4:Wab as ^ — j- 1. As soon as BB interactions are added, the starting B 
coverage is ~ 0.226 and the mixture effect disappears (Fig. 7). In the presence 
of AA and AB interactions (Fig. 8), WAB/ksT = 0.5 is sufficiently high for the 
appearance of ordered phases for A molecules, so that the resulting behavior 
is similar to that of Fig. 4. The inclusion of the three interactions (Fig. 9), 
only produces the elimination of the mixture effect, compared to Fig. 8, due 
to the low initial B coverage. 

The effect of variable BB interactions is discussed in Figs. 10 to 13. Different 
cases were considered: WAA/kBT = WAB/ksT = (Fig. 10); waa/^bT = 
and WAs/kBT = 2 (Fig. 11); WAA/kBT = 5 and WAB/ksT = (Fig. 12) and 
WAA/kBT = 5 and WAB/kBT = 2 (Fig. 13). Given the value of the parameters 
in Fig. 10, neither the coverage of A [Fig. 10 a)] nor the differential heat of ad- 
sorption [Fig. 10 c)] are affected by BB interactions. As discussed in Fig. 1, the 
coverage of B at low pressure starts at 0.5 and decreases as wsB/kBT increases 
towards the limiting value of 0.226 [see Fig. 10 a)]. However, the coverage of B 
does not present any special features. Such a special feature does indeed ap- 
pear in the behavior of qb [Fig. 10 d)], which decreases steadily as wsB/kBT 
increases below a certain critical value, w%B/kBT (being w%B/kBT ~ 2), and 
increases above it. The explanation for this behavior is that B particles adsorb 
more or less at random below WBB/kBT, thereby allowing some BB interac- 
tions which contribute to the decrease in the differential heat of adsorption. 
Above w%B/kBT, B particles adsorb forming an ordered structure so that BB 
interactions stop contributing and qb increases. The condition WAB/kBT = 
restricts the possibility of mixture effect. 

In Fig. 11, the presence of AB interactions favors the displacement of B parti- 
cles and, consequently, the slopes of the B partial isotherms are increased. On 
the other hand, as the initial fraction of B particles is high (0.3 < 6'^ < 0.5), 
which occurs for small values of WBB/kBT (0 < WBB/kBT < 0.5), AB inter- 
actions lead to mixture effect. In a wide range of coverage, 0.4 < ^ < 1, the 
curves in Fig. 11 are very similar between them (this is clearly visualized in 
the case qa and qb), which is indicative of the rapid decreasing of the number 
of B particles on the surface. 

In Fig. 12, the adsorption of A particles [Fig. 12 a)] follows a unique isotherm 
and is independent of the strength of BB interactions. Adsorption of B parti- 
cles decreases at low A coverage with increasing values oiwBB/k bT and tends 
to the limiting value 9b = 0.226. However, at high A coverage all isotherms 
tend to that corresponding to WBB/kBT = because there are no A^A^ vacant 
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sites in that range available for the adsorption of B particles. This determines 
the total coverage behavior shown in Fig. 12 b). The curves of qa-, shown in 
Fig. 12 c), present a very similar behavior between them. This is, one marked 
jump appears, corresponding to the plateau in A isotherms. In contrast, 
presents two types of behaviors [Fig. 12 d)]. Thus, at very low WBB/kBT val- 
ues (negligible interactions) and very high wsB/kBT values (above the critical 
value for the formation of the ordered phase), it remains practically constant, 
while at intermediate values it increases in line with the total coverage. In the 
last case (Fig. 13), the adsorption process can be easily understood: the B 
particles disappear for low values of iiAlksT. Then, for higher fiA/kBT's, all 
curves collapse on a unique curve. 



3.1 Phase diagrams 



In order to rationalize the results presented in previous section, we will deter- 
mine the temperature-coverage phase diagram characterizing our system in the 
range of the parameters studied. The curves will be obtained as a generaliza- 
tion of the well-known phase diagram for a lattice-gas of repulsive monomers 
adsorbed on a homogeneous square lattice, which is shown in Fig. 14. 

Some of the exact properties of this system have been found, especially by 

Onsager These are confined mostly to the special condition 6 = 1/2, but 

by symmetry, it can be deduced that 9c = 1/2, if a critical point exists, so 

this is the most interesting value of 6. Thus, the maximum of the coexistence 

curve (occurring to 9 = 1/2) corresponds to a critical value ksTc/wAA = 
1 



|21n(V2- 1) 



0.567 [29|30] . On the other hand, zones I, II and III 
correspond to a disordered lattice-gas state, a ordered state [c(2 x 2) phase], 
and a disordered lattice-hquid state, respectively. 

We start with the analysis of Fig. 13, where our model is almost identical to a 
lattice-gas of one species. As indicated in Fig. 13, the existence of repulsive AB 
interactions favors the displacement of B particles. Once 6b ~ 0, which occurs 
at iiA/kBT ^ 5, the binary mixture is equivalent to the square lattice-gas of 
one species. To corroborate this affirmation. Fig. 15 shows the total isotherms 
of Fig. 13 b) (symbols), in comparison with the adsorption isotherm corre- 
sponding to a square lattice-gas of one species and WAA/kBT = 5 (solid line). 
The analysis can be separated in two parts: for 6 < 1/2, there exist B parti- 
cles on the lattice and, consequently, the curves of the mixture deviate from 
that corresponding to one species. For 6 > 1/2, the coverage of B particles is 
negligible and all curves collapse on a unique curve. Then, the phase diagram 
characterizing a binary mixture with the set of parameters of Fig. 13 is identi- 
cal to that shown in Fig. 14. The unique difference with the one-species phase 
diagram is that the zone I now corresponds to an A-B mixture with different 
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proportions according to the value of WBs/kBT. 

As a basis for the analysis of the behavior of the system for variable wab / ^bT, 
we begin by considering one of the cases of Fig. 13 (that corresponding to 
wbe/^bT = 0), which is characterized by a phase diagram as shown in Fig. 
14. This behavior is representative of systems with high values of WAs/kBT 
(in the range WAB/kBT > 2.0), as is indicated in Fig. 16 (see the plane cor- 
responding to WAB/kBT = 2.0). As WAs/kBT is decreased, the maxima of 
the coexistence curves in zone II shift to higher values of coveragJ^. This can 
be better visualized in Fig. 17 (solid circles, bottom axis), where we plot the 
densities corresponding to the maxima of the coexistence curves in zone II, 
6c s, as a function of WAB/kBT. The figure allows to analyze the behavior of 
a binary mixture with WBs/kBT = 0, variable WABlksT and WAA/ksT in 
the critical regime. The curves can be understood according to the following 
reasoning. As it was explained for high values of the AB lateral interaction, 
the zone II corresponds to a phase of A particles. As WAB/kBT decreases, the 
phase is complemented with B particles, which are randomly distributed in 
the empty sites of the structure. As a typical example, we will analyze the 
case of voBB/ksT = WAB/kBT = 0. In this case, the B particles, which are at 
chemical potential fiB/kBT = 0, occupy at random 1/2 of the empty sites of 
the phase ^ I. Then, the phase corresponding to zone II is formed by a struc- 
ture of A particles (which occupy 1/2 of the total sites) complemented with 
B particles (which occupy at random 1/4 of the total sites). The resulting 
value of 6c is 3/4. The rest of the points in Fig. 17 can be explained by similar 
arguments. 

We now turn to the effect of BB interactions (see Fig. 17, open circles, upper 
axis). For this purpose, we set WAB/kBT = and WAA/ksT in the critical 
regime. The adsorption properties corresponding to this case were discussed 
in Fig. 12. We start the analysis with the case wsB/kBT = WAB/kBT = 0, 
where 6c = 3/4. As WBB/kBT is increased, 6c remains constant. In this case, 
the phase corresponds to a c(2 x 2) phase of A particles complemented with 
a partial coverage of B particles equal to 1/4. 

Figs. 14-17 allow to characterize the critical behavior of a binary mixture 
adsorbed in a square lattice in the region of the parameters studied in the 
present contribution. Namely: 

• In Fig. 2, the total isotherms have a pronounced plateau at 6 = 3/4 for 
strongly repulsive AA interactions, which smoothes out already for waa/ ksT < 

^ The figure shows the upper part of each coexistence curve. A complete analy- 
sis of the curve should require a more complex study (percolation of phase, zero- 
temperature calculations, etc.), which are out of the scope of the present paper. 
^ Note that 1/2 of the empty sites of the phase c(2 x 2) represents 1/4 of the total 
sites. 
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1. This result is reflected in Fig. 17. 

• Figs. 3 and 4 correspond to particular cases in Fig. 17. 

• Figs. 5 and 9 can be explained by combining the results in Fig. 17. 

• Figs. 8, 12 and 13 were discussed in details above. 

• Finally, no phase transition develops in the system when ksT/wAA > 0.567. 
This is clearly seen in Figs. 6, 7, 10 and 11, where a smooth behavior in the 
adsorption properties is obtained. 



4 Conclusions 

Using Monte Carlo simulations, we have studied the adsorption of a gas mix- 
ture of interacting particles A and B on homogeneous square surfaces. A va- 
riety of behaviors arise due to the formation of different ordered structures 
in the adlayer for different values of the lateral interactions among adsorbed 
particles. The analysis of partial and total adsorption isotherms and differen- 
tial heats of adsorption provides a detailed understanding of the adsorption 
process. This study yields to the construction of a phase-diagram which allows 
to understand the critical behavior of the system. 

An interesting feature of this work is the occurrence of a minimum in the 
global adsorption isotherm. This singularity appears as the initial fraction of 
B particles on the surface is high (0.4 < 6^^ < 0.5) and WAB/ksT > 0. This 
effect might be interesting due to the fact that seems counterintuitive to see a 
negative slope (therefore a minimum) in the total coverage. Nevertheless, the 
phenomenon is because the A particles adsorbing in the lattice expel the B 
particles at a higher ratio; then the partial A [B] coverage increases [decreases] 
(partial isotherms must cross). In the above regime, the desorbed B particles 
are more than the adsorbed A particles, therefore we see a total coverage {6) 
with a negative slope followed by a minimum. 

The computational technique used here has proven to be very powerful tool 
for these kind of lattice-gas models and many other systems in a much wider 
scope of sciences, allowing the interpretation of many experimental results 
without heavy or time-consuming calculations. 
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Figure Captions 

Fig. 1: Adsorption isotherms for the single-gas adsorption of A particles onto 
the surface showing the effect of lateral AA interactions. 

Fig. 2: Mixed-gas adsorption on a square lattice: (a) adsorption isotherms 
for A and B particles; (b) total adsorption isotherms; (c) differential heat of 
adsorption for A particles and (d) differential heat of adsorption for B particles. 
Effect of AA interactions: WAA/kBT > 0, lUBB/ksT = and WAB/ksT = 0. 

Fig. 3: As Fig. 2 for WAA/kBT > 0, wss/kBT = 2 and WAB/kBT = 0. 

Fig. 4: As Fig. 2 for WAA/ksT > 0, Wss/kBT = and WAB/kBT = 2. 

Fig. 5: As Fig. 2 for WAA/kBT > 0, WBs/kBT = 2 and WAB/kBT = 2. 

Fig. 6: Mixed-gas adsorption on a square lattice: (a) partial adsorption isotherms 
for A and B particles; (b) total adsorption isotherms; (c) differential heat of ad- 
sorption for A particles and (d) differential heat of adsorption for B particles. 
Effect of AB interactions: WAA/kBT = 0, WBB/kBT = and WAs/kBT > 0. 

Fig. 7: As Fig. 6 for WAA/kBT = 0, WBB/kBT = 2 and WAB/kBT > 0. 

Fig. 8: As Fig. 6 for WAA/ksT = 5, WBs/kBT = and WAB/kBT > 0. 

Fig. 9: As Fig. 6 for WAA/ksT = 5, WBs/kBT = 2 and WAB/kBT > 0. 

Fig. 10: Mixed-gas adsorption on a square lattice: (a) partial adsorption isotherms 
for A and B particles; (b) total adsorption isotherms; (c) differential heat of ad- 
sorption for A particles and (d) differential heat of adsorption for B particles. 
Effect of BB interactions: WAA/kBT = 0, uiBBlksT > and WAB/ksT = 0. 

Fig. 11: As Fig. 10 for WAA/kBT = 0, WBs/kBT > and WAB/kBT = 2. 

Fig. 12: As Fig. 10 for WAA/kBT = 5, WBs/kBT > and WAs/kBT = 0. 

Fig. 13: As Fig. 10 for WAA/kBT = 5, WBs/kBT > and WAs/kBT = 2. 

Fig. 14: Temperature-coverage phase diagram corresponding to a lattice-gas 
of repulsive monomers {wAA/ksT > 0) adsorbed on a homogeneous square 
lattice. 

Fig. 15: Comparison between the total isotherms in Fig. 13 b) (symbols) and 
the one corresponding to a square lattice-gas of one species with WAA/ksT = 5 
(sohd line). The inset shows a snapshot of the c(2 x 2) phase. 
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Fig. 16: Effect of the lateral interactions between A-B particles, WAB/kBT, on 
the temperature-coverage phase diagram corresponding to a binary mixture 
with WBs/kBT = and WAA/ksT in the critical regime. 

Fig. 17: Densities corresponding to the maxima of the coexistence curves in 
zone 11 as a function of (i) WAB/ksT (sohd circles, bottom axis) and (ii) wbb 
(open circles, upper axis). In case (i) [{ii)], WBs/kBT = [wAB/kBT = 0] and 
WAA/kBT is chosen in the critical regime. 
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Fig. 5d: Garcia et al. 




Fig 7c: Garcia et al. 




Fig. 14: Garcia et al. 
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Fig. 15: Garcia et al. 




Fig. 16: Garcia et al. 
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